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ABSTRACT 

The Fermi Gamma-ray Space Telescope has greatly expanded the number and energy window of observations 
of gamma-ray bursts (GRBs). However, the coarse localizations of tens to a hundred square degrees provided 
by the Fermi GRB Monitor instrument have posed a formidable obstacle to locating the bursts’ host galaxies, 
measuring their redshifts, and tracking their panchromatic afterglows. We have built a target-of-opportunity 
mode for the intermediate Palomar Transient Factory in order to perform targeted searches for Fermi afterglows. 

Here, we present the results of one year of this program: 8 afterglow discoveries out of 35 searches. Two of 
the bursts with detected afterglows (GRBs 130702A and 140606B) were at low redshift (z = 0.145 and 0.384 
respectively) and had spectroscopically confirmed broad-line Type Ic supemovae. We present our broadband 
follow-up including spectroscopy as well as X-ray, UV, optical, millimeter, and radio observations. We study 
possible selection effects in the context of the total Fermi and Swift GRB samples. We identify one new 
outlier on the Amati relation. We find that two bursts are consistent with a mildly relativistic shock breaking 
out from the progenitor star, rather than the ultra-relativistic internal shock mechanism that powers standard 
cosmological bursts. Finally, in the context of the Zwicky Transient Facility, we discuss how we will continue 
to expand this effort to find optical counterparts of binary neutron star mergers that may soon be detected by 
Advanced LIGO and Virgo. 

Subject headings: gamma-ray burst: individual (GRB 130702A, GRB 140606B) — supernovae: general — 
methods: observational — surveys — gravitational waves 
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1. INTRODUCTION 

Deep synoptic optical surveys, including the Palomar Tran¬ 
sient Factory (PTF; Law et al. 2009; Rau et al. 2009) and 
Pan-STARRS (Kaiser et al. 2010), have revealed a wealth of 
new transient and variable phenomena across a wide range 
of characteristic luminosities and timescales (Kasliwal 2011). 
With a wide (7 deg^) instantaneous field of view (FOV), mod¬ 
erately deep sensitivity (reaching R = 20.6 mag in 60 s), a con¬ 
sortium of follow-up telescopes, sophisticated image subtrac¬ 
tion and machine learning pipelines, and an international team 
of human-in-the-loop observers, PTF has been a wellspring 
of new or rare kinds of explosive transients (for instance, 
Quimby et al. 2011; Kasliwal et al. 2012) and early-time ob¬ 
servations of supernova (SNe) or their progenitors (see, for 
example, Nugent et al. 2011; Corsi et al. 2012; Ofek et al. 
2013; Gal-Yam et al. 2014). PTF has even blindly detected 
the optical emission (Cenko et al. 2014; S. B. Cenko et al., 
in preparation) from the rarest, brightest, and briefest of all 
known cosmic explosions, GRBs, hitherto only discoverable 
with the aid of precise localizations from space-based gam¬ 
ma-ray observatories. PTF has also detected explosions that 
optically resemble GRB afterglows but may entirely lack 
gamma-ray emission (Cenko et al. 2013b). 

GRBs and their broadband afterglows are notoriously chal¬ 
lenging to capture. They naturally evolve from bright to faint, 
and from high (gamma- and hard X-ray) to low (optical and 
radio) photon energies, with information encoded on energy 
scales from 1 to 10^^ GHz (Perley et al. 2014d) and timescales 
from 10“^ to 10^ s. Only with a rapid sequence of handoffs 
between facilities graded by energy passband, FOV, and po¬ 
sition accuracy have we been able to find them, pinpoint their 
host galaxies, and constrain their physics. The Swift mission 
(Gehrels et al. 2004), with its 1.4 sr-wide (50% coded) Burst 
Alert Telescope (BAT; Barthelmy et al. 2005) and its ability 
to slew and train its onboard X-ray Telescope (XRT; Burrows 
et al. 2005) and UV/Optical Telescope (UVOT; Roming et al. 
2005) on the location of a new burst within 100 s, has tri¬ 
umphed here: in nine years of operation, it has tracked down 
^700 X-ray afterglows and enabled extensive panchromatic 
observations by a worldwide collaboration of ground-based 
optical and radio facilities. 

Meanwhile, the Fermi satellite has opened up a new energy 
regime extending up to 300 GeV, with the Large Area Tele¬ 
scope (LAT; Atwood et al. 2009) detecting high-energy pho¬ 
tons for about a dozen bursts per year. The Gamma-ray Burst 
Monitor (GBM; Meegan et al. 2009), an all-sky instrument 
sensitive from 8 keV to 40 MeV, detects GRBs prolifically at 
a rate of ^ 250 yr“\ with a large number (about 44 yr“^) be¬ 
longing to the rarer short, hard bursts (Paciesas et al. 2012a). 
Although LAT can provide localizations that are as accurate 
as ^10', Fermi GBM produces error circles that are several 
degrees across. Since most bursts seem to lack GeV emission 
detectable by LAT, most Fermi GBM bursts do not receive 
deep, broadband follow-up. Consequently, their redshifts and 
the properties of their afterglows have remained largely un¬ 
known. 

As part of the intermediate Palomar Transient Factory 
(iPTF), over the past year we have developed the ability to 
rapidly tile these ~ 100 deg^ GBM error circles and pinpoint 
the afterglows. This target-of-opportunity (TOO) capability 
uses and briefiy redirects the infrastructure of our ongoing 
synoptic survey (the operation of which is discussed in Gal- 
Yam et al. 2011), notably the machine learning software and 


the instrumental pipeline composed of the Palomar 48 inch 
Oschin telescope (P48; Rahmer et al. 2008), the robotic Palo¬ 
mar 60 inch telescope (P60; Cenko et al. 2006), and associ¬ 
ated spectroscopic resources including the Palomar 200 inch 
Hale telescope (P200). 

In Singer et al. (2013b), we announced the first discov¬ 
ery of an optical afterglow based solely on a Fermi GBM 
localization.^^ That explosion, GRB 130702A/iPTF13bxl, 
was noteworthy for several reasons. First, it was detected 
by Fermi LAT. Second, it was at moderately low redshift, 
z = 0.145, yet had prompt energetics that bridged the gap be¬ 
tween “standard,” bright cosmically distant bursts and nearby 
sub-luminous bursts and X-ray hashes. Third, due to its 
low redshift, an accompanying SN was spectroscopically de¬ 
tectable. 

In this work, we begin with a detailed description of the op¬ 
eration of the iPTF GRB afterglow search. We then present 
seven more GBM-iPTF afterglows from the first 13 months 
of this project. In each of the eight cases, the association 
between the optical transient and the GRB was proven by 
the presence of high-redshift absorption lines in the opti¬ 
cal spectra and the coincident detection of a rapidly fading 
X-ray source with Swift XRT. In two cases, the positions 
were further corroborated by accurate Fermi LAT error cir¬ 
cles, and in four cases by accurate Interplanetary Network 
(IPN) triangulations involving distant spacecraft. In one case 
(GRB 140508A), the IPN triangulation was performed rapidly 
and was instrumental in selecting which optical transient can¬ 
didates to follow up. In six cases, radio afterglows were de¬ 
tected. Our discovery rate of 8 out of 35 events is consistent 
with the ages and searched areas of the GBM bursts, com¬ 
bined with the luminosity function of optical afterglows. Con¬ 
sequently, by tiling larger areas and/or stacking exposures, the 
iPTF afterglow search should be able to scale to coarser local¬ 
izations and fainter optical signals, such as those associated 
with short GRBs. 

Next, we present extensive follow-up observations, includ¬ 
ing R-band photometry from the P48, multicolor photometry 
from the P60, spectroscopy (acquired with the P200, Keck, 
Gemini, APO, Magellan, Very Large Telescope (VLT), and 
GTC), and radio observations with the Karl G. Jansky Very 
Large Array(VLA), the Combined Array for Research in 
Millimeter-wave Astronomy (CARMA; Bock et al. 2006; 
Corder et al. 2010), the Australia Telescope Compact Array 
(ATCA; Frater et al. 1992), and the Arcminute Microkelvin 
Imager (AMI; Zwart et al. 2008). We provide basic physi¬ 
cal interpretations of the broadband spectral energy distribu¬ 
tions (SEDs) of these afterglows. We find that seven of the 
events are consistent with the classic model of synchrotron 
cooling of electrons that have been accelerated by a single 
forward shock encountering either the constant-density cir- 
cumburst interstellar medium (ISM; broadband behavior pre¬ 
dicted in Sari et al. 1998) or a stellar (i.e., Wolf-Ray et) wind 
environment (Chevalier & Li 1999). The possible excep¬ 
tion, GRB 140620A/iPTF14cva, can probably be explained 
by standard extensions of this model, a reverse shock or an 
inverse Compton component. 

There are two earlier related cases. The optical afterglow of 
GRB 090902B was detected ex post facto in tiled observations with Robotic 
Optical Transient Search (ROTSE) about 80 minutes after the burst, but the 
afterglow was initially discovered with the help of an X-ray detection in Swift 
observations of the LAT error circle. GRB 120716A was identified by iPTF 
by searching a 2 deg^ IPN error box (Cenko et al. 2012). 
http://www.vla.nrao.edu 
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Two of the afterglows (GRB 130702A/iPTF13bxl and 
GRB 140606B/iPTF14bfu) faded away to reveal spectroscop¬ 
ically detected broad-line Type Ic SNe (SNelc-BL). Despite 
the abundant photometric evidence for SNe in afterglow light 
curves (see Li & Hjorth 2014 and references therein), the 
distinction of SN spectroscopy has been shared by scarcely 
tens^^ out of ^800 long Swift bursts in nine years of opera¬ 
tion. 

We estimate the kinetic energies of the relativistic blast 
waves of these events from their X-ray afterglows (Freed¬ 
man & Waxman 2001). We find that although the 
gamma-ray energetics of these eight bursts are broadly 
similar to the Swift sample, two low-luminosity bursts 
(GRBs 130702A and 140606B) have significantly lower ki¬ 
netic energies. We discuss the possibility that these two bursts 
arise not from a standard ultra-relativistic internal shock, but 
from a mildly relativistic shock as it breaks out from the pro¬ 
genitor star (see, for example, Nakar & Sari 2012). 

We conclude by discussing prospects for targeted optical 
transient searches in wide areas. This is especially relevant 
for optical counterparts of gravitational wave (GW) events. 
We illustrate that optical afterglows of short bursts, which are 
intimately linked to the prime sources for the Advanced Laser 
Interferometer GW Observatory (LIGO) and Virgo, should 
be well within the reach of a similar approach using Zwicky 
Transient Facility (ZTF; Kulkami 2012; Bellm 2014; Smith 
et al. 2014). 


2. SEARCH METHODOLOGY 

We begin by describing our TOO observations and after¬ 
glow search step by step. 

2.1. Automated TOO Marshal: Alerts and Tiling 

A program called the iPTE TOO Marshal moni¬ 
tors the stream of Gamma-ray Coordinates Network 
(GCN) notices^^ from the three redundant, anonymous 
NASA/GSEC VOEvent servers. It listens for notices of 
type FERMI_GBM_GND_POS, sent by GBM’s automated 
on-ground localization, or FERMI_GBM_FIN_POS, sent by 
the GBM burst advocate’s human-in-the-loop localization.^^ 

Upon receiving either kind of notice, the TOO Marshal de¬ 
termines if the best-estimate sky position is observable from 
Palomar at any time within the 24 hr after the trigger. The 
criterion for observability is that the position is at an altitude 
> 23.°5 (i.e. airmass < 2.5), at least 20° from the center of 
the moon, at an hour angle between ±6!^5, and that the Sun is 
at least 12° below the horizon at Palomar. 

If the position is observable and the Icr statistical error ra¬ 
dius Tstat reported in the GCN notice is less than 10°, the 
TOO Marshal selects a set of 10 P48 fields that optimally 
cover the error region.It converts the GBM position es¬ 
timate and radius into a probability distribution by applying 

Between photometric, late-time red bumps and unambiguous spectral 
identifications, there are also GRB-SNe that have some SN-associated spec¬ 
tral features. The number of GRBs with spectroscopic SNe is, therefore, ill 
defined. See Hjorth & Bloom (2012, p. 169, and references therein) for a 
more complete census. 

http://gcn.gsfc.nasa.gov 

Usually, the Fermi team suppresses the notices if the burst is detected 
and localized more accurately by Swift BAT. 

We made one exception to our GBM error radius cutoff: we followed 
up GRB 140219A, which had a GBM error circle with a radius of 12.8°, but 
had an IPN localization spanning 0.6 deg^ (Hurley et al. 2014a). Despite 
searching about 80% of the IPN polygon, we detected no afterglow (Singer 
et al. 2014d). This is potentially a dark burst candidate. 


a well-known empirical prescription of the systematic errors 
of the GBM localization. Paciesas et al. (2012b) state that 
the total effective error radius in the FERMI_GBM_FIN_POS 
localizations is well described by the quadrature sum of the 
statistical radius and a systematic contribution, where the 
systematic is 2.°6 for 72% of bursts and 10.°4 for 28% of 
bursts. We use the weighted rms of these two values, r^y^ = 

Y^0.72(2.°6)^-1-0.28(10.°4)2 « 6°. The total error radius is 

then Teff = \/^stat^ + ^sys^- Wc construct a Eisher-von Mises 
distribution, centered on the best-estimate position, with a 
concentration parameter of 


K, = 


1 - COS 



( 1 ) 


With the FERMI_GBM_FIN_POS alert, the Fermi GBM 
team also distributes a detailed localization map that accounts 
for the systematic effects (Connaughton et al. 2015). The 
TOO Marshal retrieves from the Fermi data archive a file 
that describes the la, 2a, and 3a significance contours. If 
the localization has significant asymmetry, we also retrieve a 
2D FITS image whose pixel values correspond to the GBM 
localization significance, and use this instead of the Fish- 
er-von Mises distribution. 

Giving preference to fields for which deep co-added ref¬ 
erence images exist, the TOO Marshal selects 10 P48 fields 
spanning an area of ^ 72 deg^ to maximize the probability of 
enclosing the true (but as yet unknown) location of the source, 
assuming the above distribution. 

The Marshal then immediately contacts a team of humans 
(the authors) by SMS text message, telephone, and e-mail. 
The humans are directed to a mobile-optimized web applica¬ 
tion to trigger the P48 (see Fig. 11). 


2.2. Triggering the P48 

Within the above constraints, we decide whether to follow 
up the burst based on the following criteria. The event must 
be <12 hr old when it first becomes observable from Palomar, 
and we must cover enough of the error circle to have a ^30% 
chance of enclosing the position of the source. We discard 
any bursts that are detected and accurately localized by Swift 
BAT, because these are more efficiently followed up by con¬ 
ventional means. We also give preference to events that are 
out of the Galactic plane and that are observable for at least 
3 hr. 

There are some exceptional circumstances that override 
these considerations. If the burst’s position estimate is acces¬ 
sible within an hour after the burst, we may select it even if 
the observability window is very brief. If the burst is very well 
localized or has the possibility of a substantially improved lo¬ 
calization later due to a LAT or IPN detection, we may select 
it even if it is in the Galactic plane. 

The default observing program is three epochs of P48 im¬ 
ages at a 30-minute cadence. The human may shorten or 
lengthen the cadence if the burst is very young or old (see 
the discussion of Equation (2) in Section 2.4 below), change 
the number of epochs, or add and remove P48 fields. When 
the human presses the “Go” button, the TOO Marshal sends a 
machine-readable e-mail to the P48 robot. The robot adds the 
requested fields to the night’s schedule with the highest pos¬ 
sible priority, ensuring that they are observed as soon as they 
are visible. 


2.3. Automated Candidate Selection 
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As the night progresses, the TOO Marshal monitors the 
progress of the observations and the iPTF real-time image 
subtraction pipeline (R E. Nugent et al. 2015, in preparation). 
The real-time pipeline creates difference images between the 
new P48 observations and co-added references composed of 
observations from months or years earlier. It generates can¬ 
didates by performing source extraction on the difference im¬ 
ages. A machine learning classifier assigns a real/bogus score 
(RB2; Brink et al. 2013) to each candidate that predicts how 
likely the candidate is to be a genuine astrophysical source 
(rather than a radiation hit, a ghost, an imperfect image sub¬ 
traction residual, or any other kind of artifact). 

Table 1 lists the number of candidates that remain after 
each stage of candidate selection. First, requiring candidates 
to have signal-to-noise ratio (S/N)> 5 gives us a median of 
35,000 candidates. This number varies widely with galac¬ 
tic latitude and the area searched (a median of ^500 deg“^). 
Second, we only select candidates that have RB2> 0.1, re¬ 
ducing the number of candidates to a median of 36% of the 
original list.'^^ Third, we reject candidates that coincide with 
known stars in reference catalogs (Sloan Digital Sky Sur¬ 
vey (SDSS) and the PTE reference catalog), cutting the list 
to 17%. Fourth, we eliminate asteroids cataloged by the Mi¬ 
nor Planet Center, reducing the list to 16%. Fifth, we demand 
at least two secure P48 detections after the GBM trigger, re¬ 
ducing the list to a few percent, or ^ 500 candidates. 

When the image subtraction pipeline has finished analyzing 
at least two successive epochs of any one field, the TOO Mar¬ 
shal contacts the humans again and the surviving candidates 
are presented to the humans via the Treasures portal. 

2.4. Visual scanning in Treasures Portal 

The remaining candidate vetting steps currently involve 
human participation and are informed by the nature of the 
other transients that iPTF commonly detects: foreground SNe 
(slowly varying and in low-z host galaxies), active galactic 
nuclei (AGNs), cataclysmic variables, and M-dwarf flares. 

In the Treasures portal, we visually scan through the auto¬ 
matically selected candidates one P48 field at a time, exam¬ 
ining ^10 objects per field (see Figure 12 for a screenshot of 
the Treasures portal). We visually assess each candidate’s im¬ 
age subtraction residual compared to the neighboring stars of 
similar brightness in the new image. If the residual resembles 
the new image’s point-spread function, then the candidate is 
considered likely to be a genuine transient or variable source. 

Next, we look at the photometric history of the candidates. 
Given the time, t, of the optical observation relative to the 
burst and the cadence, dt, we expect that a typical optical 
afterglow that decays as a power law Tv with a = 1, 

would fade by 5m = 2.5\ogiQ{\+5t/t) mag over the course of 
our observations. Any source that exhibits statistically signif¬ 
icant fading {5m jm ^1) consistent with an afterglow decay 
becomes a prime target."^^ 

Note that ala decay in brightness requires such a source to 
be 

-2.51ogio(^) (2) 

This RB2 threshold is somewhat deeper than that which is used in the 
iPTF survey. An improved classifier, RB4 (Bue et al. 2014), entered evalua¬ 
tion in 2014 August shortly before GRB 140808A. 

A source that exhibits a statistically significant rise is generally also 
followed up, but as part of the main iPTF transient survey, rather than as a 
potential optical afterglow. 


brighter than the la limiting magnitude of the exposures. 
For example, given the P48’s typical limiting magnitude of 
R = 20.6 and the standard cadence of = 0.5 hr, if a burst is 
observed f = 3 hr after the trigger, its afterglow may be ex¬ 
pected to have detectable photometric evolution only if it is 
brighter than R = 18.3. Noting that long GRBs preferentially 
occur at high redshifts and in intrinsically small, faint galaxies 
(Svensson et al. 2010), we consider faint sources that do not 
display evidence of fading if they are not spatially coincident 
with any sources in SDSS or archival iPTF observations. 

If a faint source is near a spatially resolved galaxy, then 
we compute its distance modulus using the galaxy’s redshift 
or photometric redshift from SDSS. We know that long GRB 
optical afterglows at f = 1 day typically have absolute mag¬ 
nitudes of -25 < Mb < -21 mag (la range; see Figure 9 
of Kann et al. 2011). Most SNe are significantly fainter: 
Type la are typically Mb ^ -19 mag whereas Ibc and II are 
Mb ^ -17 mag, with luminous varieties of both Type Ibc 
and II extending to Mb ^ -19 mag (Richardson et al. 2002; 
Li et al. 2011). Therefore, if the candidate’s presumed host 
galaxy would give it an absolute magnitude Mr < -20 mag, it 
is considered promising. This criterion is only useful for long 
GRBs because short GRB afterglows are typically ^ 6 mag 
fainter than long GRB afterglows (Kann et al. 201 1). 

The human saves all candidates that are considered promis¬ 
ing by these measures to the iPTF Transient Marshal database. 
This step baptizes them with an iPTF transient name, which 
consists of the last two digits of the year and a sequential al¬ 
phabetic designation. 

2.5. Archival vetting in the Transient Marshal 

Once named in the Transient Marshal, we perform archival 
vetting of each candidate using databases including VizieR 
(Ochsenbein et al. 2000), NED,"^^ the High Energy Astro¬ 
physics Science Archive Research Center (HEASARC),'^^ 
and Catalina Real-time Transient Survey (Drake et al. 2009), 
in order to check for any past history of variability at that posi¬ 
tion (see Figure 13 for a screenshot of the Transient Marshal). 

We check for associations with known quasars or AGNs 
in Veron-Cetty & Veron (2010) or with AGN candidates in 
Flesch (2010). 

M dwarfs can produce bright, blue, rapidly fading opti¬ 
cal flares than can mimic optical afterglows. To filter our 
M dwarfs, we check for quiescent infrared counterparts in 
WISE (Cutri & et al. 2014). Stars of spectral type L9-M0 
peak slightly blueward of the WISE bandpass, with typical 
colors (Wright et al. 2010) 

3< [R-Wl] <12 
0.1 <[W1-W2]<0.6 
0.2<[W2-W3]<1 
0<[W3-W4]<0.2. 

Therefore, a source that is detectable in WISE but that is ei¬ 
ther absent from or very faint in the iPTF reference images 
suggests a quiescent dwarf star. 

2.6. Photometric, Spectroscopic, and Broad-band 
Eollow-up 

The above stages usually result in ^10 promising optical 
transient candidates that merit further follow-up. If, by this 

http://ned.ipac.caltech.edu 

http://heasarc.gsfc.nasa.gov 
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Table 1 

Number of Optical Transient Candidates Surviving Each Vetting Stage 


GRB 

S/N 

>5 

RB2 

>0.1 

Not 

Stellar 

Not in 

MPC^ 

Detected 

Twice 

Saved for 
Follow-up 

RB2'’ 

130702A 

14629 

2 388 

1346 

1323 

417 

11 

0.843 

131011A 

21308 

8 652 

4344 

4197 

434 

23 

0.198 

131231A 

9 843 

2503 

1776 

1543 

1265 

10 

0.137 

140508A 

48747 

22673 

9970 

9969 

619 

42 

0.730 

140606B 

68 628 

26070 

11063 

11063 

1449 

28 

0.804 

140620A 

152224 

50930 

17 872 

17 872 

1904 

34 

0.826 

140623A 

71219 

29434 

26279 

26279 

442 

23 

0.873 

140808A 

19 853 

4 804 

2 349 

2 349 

79 

12 

0.318 


Median reduction 

36% 

17% 

16% 

1.7% 

0.068% 



^ Not in Minor Planet Center database 
^ RB2 score of optical afterglow in earliest P48 detection 


point, data from Fermi LAT or from IPN satellites are avail¬ 
able, we can use the improved localization to select an even 
smaller number of follow-up targets. 

For sources whose photometric evolution is not clear, we 
perform photometric follow-up. We may schedule additional 
observations of some of the P48 fields if a significant num¬ 
ber of candidates are in the same field. We may also use 
the P48 to gather more photometry for sources that are su¬ 
perimposed on a quiescent source or galaxy, in order to make 
use of the image subtraction pipeline to automatically obtain 
host-subtracted magnitudes. For isolated sources, we sched¬ 
ule one or more epochs of r-band photometry with the P60. If, 
by this point, any candidates show strong evidence of fading, 
we begin multicolor photometric monitoring with the P60. 

Next, we acquire spectra for one to three candidates per 
burst using the P200, Gemini, Keck, Magellan, or Himalayan 
Chandra Telescope (HCT). A spectrum that has a relatively 
featureless continuum and high-redshift absorption lines se¬ 
cures the classification of the candidate as an optical after¬ 
glow. 

Once any single candidate becomes strongly favored over 
the others based on photometry or spectroscopy, we trigger 
X-ray and UV observations with Swift and radio observations 
with VLA, CARMA, and AMI. Detection of a radio or X-ray 
afterglow typically confirms the nature of the optical transient, 
even without spectroscopy. 

Finally, we promptly release our candidates, upper limits, 
and/or confirmed afterglow discovery in GCN circulars. 

2.7. Long-term Monitoring and Data Reduction 

The reported P48 magnitudes are all in the Mould R band 
and in the AB system (Oke & Gunn 1983), calibrated with 
respect to either r' point sources from SDSS or for non-SDSS 
fields using the methods described in Ofek et al. (2012). 

To monitor the optical evolution of afterglows identified by 
our program, we typically request nightly observations in ri 
(and occasionally gz) filters for as long as the afterglow re¬ 
mained detectable. Bias subtraction, fiat-fielding, and other 
basic reductions are performed automatically at Palomar by 
the P60 automated pipeline using standard techniques. Im¬ 
ages are then downloaded and stacked as necessary to im¬ 
prove the S/N. Photometry of the optical afterglow is then 
performed in IDL using a custom aperture-photometry rou¬ 
tine, calibrated relative to SDSS secondary standards in the 
field (when available) or using our own solution for secondary 
field standards constructed during a photometric night (for 
fields outside the SDSS footprint). 


For some bursts (GRB 140606B), we also obtain pho¬ 
tometry with the Large Monolithic Imager (LMI) mounted 
on the 4.3 m Discovery Channel Telescope (DCT) in Happy 
Jack, Arizona. Standard CCD reduction techniques (e.g., 
bias subtraction, fiat-fielding) are applied using a custom 
IRAF pipeline. Individual exposures are aligned with re¬ 
spect to astrometry from the Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006) using SCAMP (Bertin 2006) 
and stacked with SWarp (Bertin et al. 2002). 

Where GROND (Greiner et al. 2008) and RATIR (Butler 
et al. 2012) have reported multicolor photometry in GCN cir¬ 
culars, we include their published data in Table 5 and our 
light-curve plots. 

We monitor GBM-iPTF afterglows with CARMA, a mil¬ 
limeter-wave interferometer located at Cedar Flat near Big 
Pine, California. All observations are conducted at 93 GHz 
in single-polarization mode in the array’s C, D, or E configu¬ 
ration. Targets are typically observed once for 1-3 hr within a 
few days after the GRB, establishing the phase calibration us¬ 
ing periodic observations of a nearby phase calibrator and the 
bandpass and the fiux calibration by observations of a stan¬ 
dard source at the start of the track. If detected, we acquire ad¬ 
ditional observations in approximately logarithmically spaced 
time intervals until the afterglow fiux falls below detection 
limits. All observations are reduced using MIRIAD using 
standard interferometric flagging and cleaning procedures. 

We look for radio afterglows at 6.1 and/or 22 GHz with 
VLA. VLA observations are reduced using the Common As¬ 
tronomy Software Applications (CASA) package. The cali¬ 
bration is performed using the VLA calibration pipeline. Af¬ 
ter running the pipeline, we inspect the data (calibrators and 
target source) and apply further fiagging when needed. The 
VLA measurement errors are a combination of the rms map 
error, which measures the contribution of small unresolved 
fiuctuations in the background emission and random map fiuc- 
tuations due to receiver noise, and a basic fractional error 
(here estimated to be ^5%) which accounts for inaccura¬ 
cies of the fiux density calibration. These errors are added 
in quadrature, and total errors are reported in Table 6. 

Starting in 2014 August, we also look for radio emission 
with AMI. AMI is composed of eight 12.8 m dishes operating 
in the 13.9-17.5 GHz range (central frequency of 15.7 GHz) 
when using frequency channels 3-7 (channels 1, 2, and 8 are 
disregarded due to their currently susceptibility to radio in¬ 
terference). For further details on the reduction and analysis 
performed on the AMI observations please see Anderson et al. 
(2014b). 
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3. THE GBM-iPTF BURSTS 

To date, we have successfully followed up 35 Fermi GBM 
bursts and detected eight optical afterglows. The detections 
are listed in Table 2, and all of the P48 tilings are listed in Ta¬ 
ble 3. Figure 5 shows the GBM localizations and P48 tilings 
for the detected bursts. In Figure 1, the light curves are shown 
in the context of a comprehensive sample of long GRB af¬ 
terglows compiled by D. A. Kann (2015, private communica¬ 
tion). 

The outcome of an individual afterglow search is largely de¬ 
termined by two factors: how much probability is contained 
within the P48 footprints, and how bright the afterglow is at 
the time of the observations (see Figure 2). We calculate the 
expected success rate as follows. For each burst, we find the 
prior probability that the position is contained within the P48 
fields that we observed. We then compute the fraction of af¬ 
terglows from Kann’s sample (which has a mean and standard 
deviation of 22 ± 2 mag at ^ = 1 day) that are brighter than 
R = 20.6 mag at the same age as when the P48 observations 
started. The product of these two numbers is the prior proba¬ 
bility of detection for that burst. By summing over all of the 
iPTF/GBM bursts, we obtain the expected number of detec¬ 
tions. Within 95% confidence bootstrap error bars, we find an 
expected 5.5-8.5 detections, or a success rate of 16%-24%. 
This is consistent with the actual success rate of 23%. 

This suggests that the success rate is currently limited by 
the survey area and the response time (dictated by sky posi¬ 
tion and weather). We could increase the success rate by de¬ 
creasing the maximum time since trigger at which we begin 
follow-up. We could increase the success rate without ad¬ 
versely affecting the number of detections by simply search¬ 
ing a greater area for coarsely localized events. 

Fig. Set 3. Light curves and SEDs 

Over the next few sections, we summarize the observations 
and general physical interpretation of all of the GBM-iPTF 
afterglows detected to date. Figure 3 shows the light curves 
and SEDs spanning X-ray, UV, optical, IR, and radio frequen¬ 
cies. Table 4 contains a log of our spectroscopic observations. 
Table 5 lists a selection of ultraviolet, optical, and infrared 
observations, including all of our P48 and P60 observations. 
Table 6 lists all of our radio detections. 

3.1. GRB 130702A/iPTF13bxl 

This is the first GBM burst whose afterglow we discovered 
with iPTF (Singer et al. 2013b), indeed the first afterglow 
ever to be pinpointed based solely on a Fermi GBM local¬ 
ization. It is also the lowest-redshift GRB in our sample, so 
it has the richest and most densely sampled broadband after¬ 
glow data. It has two other major distinctions: its associated 
SN (SN 2013dx, Schulze et al. 2013; Pozanenko et al. 2013; 
Cenko et al. 2013a; D’Elia et al. 2013) was detected spectro¬ 
scopically, and its prompt energetics are intermediate between 
low-luminosity GRBs (llGRBs) and standard cosmic bursts 
(see below). 

Based on the Fermi GBM ground localization with an er¬ 
ror radius of 4°, we imaged 10 fields twice with the P48 
at Ar = ^-^gbm = 4.2 hr after the burst.^"^ We scheduled 

At the time, our tiling algorithm selected fields based on an empirical 
calibration of Fermi GBM’s systematic errors. We had selected bursts that 
were detected by both Swift and Fermi and constructed a fit to a cumulative 
histogram of the number of bursts whose BAT or XRT positions were within 
a given number of nominal Icr statistical radii of the center of the Fermi error 
circle. Our tiling algorithm scaled this fit by the la radius of the burst in 


P60 imaging and P200 spectroscopy for three significantly 
varying sources. Of the three, iPTF13bxl showed the clear¬ 
est evidence of fading in the P48 images. Its spectrum at 
Af = 1.2 days consisted of a featureless blue continuum. We 
triggered Swift, which found a bright X-ray source at the posi¬ 
tion of iPTFlSbxl (Singer et al. 2013c; D’Avanzo et al. 2013). 
Shortly after we issued our GCN circular (Singer et al. 2013c), 
Cheung et al. (2013) announced that the burst had entered the 
FOV of FAT at Af = 250 s. The FAT error circle had a radius 
of 0.°5, and its center was 0.°8 from iPTF13bxl. An IPN trian¬ 
gulation Wi\h MESSENGER (GRNS), INTEGRAL (SPI-ACS), 
E^rm/-GBM, and Konus-WmJ (Hurley et al. 2013) yielded a 
0.°46-wide annulus that was also consistent with the OT. 

The afterglow’s position is Of'6 from an R = 23.01 mag 
source that is just barely discernible in the P48 reference im¬ 
ages. A spectrum from NOT-f ALFOSC (Leloudas et al. 2013) 
determined a redshift of z = 0.145 for a galaxy 7."6 to the 
south of iPTF13bxl. At Af = 2.0 days, we obtained a Magel- 
lan-flMACS spectrum (Mulchaey et al. 2013) and found weak 
emission lines at the location of the afterglow that we inter¬ 
preted as Ha and [O III] at the same redshift. Kelly et al. 
(2013) characterized the burst’s host environment in detail and 
concluded that it exploded in a dwarf satellite galaxy. 

Joining the two P48 observations at Ar < 1 day to the 
late-time P60 light curve requires a break at At = 1.17 ± 
0.09 days, with slopes (To,i = 0.57 ±0.03 and (To ,2 = 1.05 ± 
0.03 before and after the break, respectively. The XRT 
light-curve begins just prior to this apparent break and seems 
to follow the late-time optical decay (until the SN begins to 
dominate at Ar = 5 days), although the automated Swift light 
curve analysis (Evans et al. 2009) also suggests a possible 
X-ray break with about the same time and slopes. This hints 
at an achromatic break, normally a signature of a jet. How¬ 
ever, the late slope and the change in slope are both unusually 
shallow for a jet break. Furthermore, the radio light curve 
does not exhibit a break. The change in slope is also a little 
too large for cooling frequency crossing the band (for which 
one would expect Aa =1/4). An energy injection or a struc¬ 
tured jet model may provide a better fit (Panaitescu 2005). 

Late-time At > I day observations include several P60 gri 
observations, three RATIR r'i'ZYJH epochs, an extensive 
Swift XRT and UVOT light curve, and radio observations with 
VLA and CARMA (although of the VLA data, we only have 
access to the first observation). The optical and X-ray spec¬ 
tral slopes are similar, /3o = 0.7 ±0.1 and px = 0.8 ± 0.1. An 
SED at 2 < Ar < 2.3 days is well explained by the standard 
external shock model (Sari et al. 1998) in the slow cooling 
regime, with lying between the VLA and CARMA fre¬ 
quencies and Uc in the optical. This fit requires a relatively 
fiat electron spectrum, dn^/d^Q oc 7e“^ with p 1.6, cut off 
at high energies. Applying the relevant closure relations (for 
the case of 1 < p < 2, see Dai & Cheng 2001) to ax and Px 
permits either an ISM or wind environment. 

Our late-time spectroscopy and analysis of the SN will be 
published separately (S. B. Cenko et al. 2015, in preparation). 

3.2. GRB 13101 lA/iPTF13dsw 


question and then constructed a 2D angular probability distribution from it. 
For sufficiently large error radii, this prescription produced probability dis¬ 
tributions that had a hole in the middle. For this reason, the tiling algorithm 
picked out P48 fields that formed an annulus around the GBM Icr error cir¬ 
cle (not, as we stated in Singer et al. 2013b, because of a lack of reference 
images). 
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Figure 1. Optical light curves of Fermi-\PT¥ afterglows to date. The light curves of the eight iPTF/GBM bursts are shown in red. For comparison, the gray lines 
show a comprehensive sample of long GRB optical light curves from Cenko et al. (2009), Kann et al. (2010), Perley et al. (2014d), and D. A. Kann (2015, private 
communication). The white area outside of the light-gray shading illustrates the range of GRB afterglows that are accessible given a half-hour cadence and the 
P48’s 60 s limiting magnitude of R = 20.6. The two light curves shown in blue are other related iPTF transients. The hrst is PTFl lagg, an afterglow-like transient 
with no detected gamma-ray emission (Cenko et al. 2013b). The second is GRB 140226A/iPTF14yb, reported initially by iPTF from its optical afterglow (Cenko 
et al. 2014; Cenko et al., in preparation), and later by IPN from its gamma-ray emission (Hurley et al. 2014b). 


We started P48 observations of Fermi trigger 403206457 
(Jenke 2013) about 11.6 hr after the burst. The optical tran¬ 
sient iPTF13dsw (Kasliwal et al. 2013) faded from R = 
19.7 mag io R = 20.2 mag from 11.6 to 14.3 hr. The latest 
pre-trigger image on 2013 September 25 had no source at this 
location to a limit of /? > 20.6 mag. The optical transient con¬ 
tinued to fade as it was monitored by several facilities (Xu 
et al. 2013a, c; Perley et al. 2013; Sudilovsky et al. 2013; Voi¬ 
nova et al. 2013). 


At 15.1 hr after the burst, we obtained a spectrum 
of iPTF13dsw with the Gemini Multi-object Spectrograph 
(GMOS) on the Gemini-South telescope. GMOS was con¬ 
figured with the R400 grating with a central wavelength of 
7200 A and the F'slit, providing coverage over the wave¬ 
length range of 5100-9300 A with a resolution of 3 A. No 
prominent features were detected over this bandpass, while 
the spectrum had a typical SNR of ^ 3-4 per 1.4 A pixel. Rau 
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Table 2 

GBM-iPTF Detections 


GRB 

OT 

R.A. 

(J2000) 

Deck 

(J2000) 

Gal. 

Lat.^ 

z 

•^peak 

(keV, rest) 

■^ 7 , iso 

(10^^ erg, rest)^’° 

^90 (s) 


GRB 130702A 

iPTF13bxl 

14'’29"'15® 

+15°46'26" 

65° 

0.145 

18± 3 

<0.065± 0.001 

58.9± 

6.2 

17.38 

GRB 131011A 

iPTF13dsw 

02'’10”06* 

.4024'40" 

-61° 

1.874 

625± 92 

14.606± 1.256 

77.1± 

3 

19.83 

GRB 131231A 

iPTF13ekl 

00'’42”22" 

-1°39'11" 

-64° 

0.6419 

291± 6 

23.015± 0.278 

31.2± 

0.6 

15.85 

GRB 140508A 

iPTF14aue 

17'’01”52" 

+46°46'50" 

38° 

1.03 

534± 28 

24.529± 0.86 

44.3± 

0.2 

17.89 

GRB 140606B 

iPTF14bfu 

2l'’52”30* 

+32°00'51" 

-17° 

0.384 

801± 182 

0.468± 0.04 

22.8± 

2.1 

19.89 

GRB 140620A 

iPTF14cva 

18'’47"'29" 

+49043/52" 

21° 

2.04 

387± 34 

7.28 ± 0.372 

45.8± 

12.1 

17.60 

GRB 140623A 

iPTF14cyb 

15'’01"'53® 

+81°11'29" 

34° 

1.92 

834± 317 

3.58 ± 0.398 

114.7± 

9.2 

18.04 

GRB 140808A 

iPTF14eag 

14'’44"'53® 

+49°12'51" 

59° 

3.29 

503± 35 

8.714± 0.596 

4.5± 

0.4 

19.01 


^ Galactic latitude of optical afterglow. This is one of the main factors that influences the number of optical transient candidates in Table 1 . 

^ is given for a 1 keV-lOMeV rest-frame bandpass. 

^ The rest-frame spectral properties, £’peak and for GRB 130702A are reproduced from Amati et al. (2013). For all other bursts, we calculated these 

quantities from the spectral fits (the scat files) in the Fermi GBM catalog (Goldstein et al. 2012) using the k-correction procedure described by Bloom et al. 

( 2001 ). 

^ R-band apparent magnitude in initial P48 detection. 


Table 3 

Log of P48 Tilings for Fermi GBM Bursts 


GRB Time^ 

GBM 

Fluence*^ 

^P48 

“^burst 

P48 

Area^ 

Prob.° 

2013 Jun 28 20:37:57 

10 ±0.1 

10.02 

73 

32% 

^2013 Jul 02 00:05:20 

57 ±1.2 

4.20 

74 

38% 

2013 Aug 28 07:19:56 

372 ±0.6 

20.28 

74 

64% 

2013 Sep 24 06:06:45 

37 ±0.6 

23.24 

74 

28% 

2013 Oct 06 20:09:48 

18 ±0.6 

15.26 

74 

18% 

^2013 Oct 11 17:47:30 

89 ±0.6 

11.56 

73 

54% 

2013 Nov 08 00:34:39 

28 ±0.5 

4.69 

73 

37% 

2013 Nov 10 08:56:58 

33 ±0.3 

17.47 

73 

44% 

2013 Nov 25 16:32:47 

5.5±0.3 

11.72 

95 

26% 

2013 Nov 26 03:54:06 

17 ±0.3 

6.94 

109 

59% 

2013 Nov 27 14:12:14 

385 ±1.4 

13.46 

60 

50% 

2013 Dec 30 19:24:06 

41 ±0.4 

7.22 

80 

38% 

^2013 Dec 31 04:45:12 

1519 ±1.2 

1.37 

30 

32% 

2014 Jan 04 17:32:00 

333 ±0.6 

18.57 

15 

11% 

2014 Jan 05 01:32:57 

6.4±0.1 

7.63 

74 

22% 

2014 Jan 22 14:19:44 

9.1±0.5 

11.97 

75 

34% 

2014 Feb 11 02:10:41 

7.4±0.3 

1.77 

44 

19% 

2014 Feb 19 19:46:32 

28 ±0.5 

7.01 

71 

14% 

2014 Feb 24 18:55:20 

24 ±0.6 

7.90 

72 

30% 

2014 Mar 11 14:49:13 

40 ±1.2 

12.18 

73 

54% 

2014 Mar 19 23:08:30 

71 ±0.3 

3.88 

74 

48% 

2014 Apr 04 04:06:48 

82 ±0.2 

0.11 

109 

69% 

2014 Apr 29 23:24:42 

6.2±0.2 

10.99 

74 

15% 

^2014 May 08 03:03:55 

614 ±1.2 

6.68 

73 

67% 

2014 May 17 19:31:18 

45 ±0.4 

8.60 

95 

69% 

2014 May 19 01:01:45 

39 ±0.5 

4.42 

73 

41% 

^2014 Jun 06 03:11:52 

76 ±0.4 

4.08 

74 

56% 

2014 Jun 08 17:07:11 

19 ±0.6 

11.20 

73 

49% 

^2014 Jun 20 05:15:28 

61 ±0.6 

0.17 

147 

59% 

^2014 Jun 23 05:22:07 

61 ±0.6 

0.18 

74 

4% 

2014 Jun 28 16:53:19 

18 ±1.0 

16.16 

76 

20% 

2014 Jul 16 07:20:13 

2.4±0.3 

0.17 

74 

28% 

2014 Jul 29 00:36:54 

81 ±0.7 

3.43 

73 

65% 

2014 Aug 07 11:59:33 

13 ±0.1 

15.88 

73 

54% 

^2014 Aug 08 00:54:01 

32 ±0.3 

3.25 

95 

69% 


^ Time of Fermi GBM trigger. —)■ Afterglow detections are marked with an 
arrow. The corresponding entries in Table 2 can be found by matching the 
date to the GRB name (GRB YYMMDDA). 

^ Observed Fermi GBM fluence in the 10-1000 keV band, in units 
of 10“^ erg cm“^. This quantity is taken from the beat files from the Fermi 
GRB catalog at HE AS ARC. 

^ Age in hours of the burst at the beginning of the P48 observations. 

^ Area in deg^ spanned by the P48 fields. 

^ Probability, given the Fermi GBM localization, that the source is contained 
within the P48 fields. 


1.0 

-^^^- 

1- 

0.8 

- 

- 


• 14eag^j4aue 


0.6 

|14cva r- 


0.4 

#13bxl 

- 


#13ekl 


0.2 

- 

- 

0.0 

|14cyb 1 1 1 

J_ 


age at first obervation (days) 


Figure 2. Prior probability of containing the burst’s location within the P48 
fields vs. age of the burst at the beginning of P48 observations. Afterglow 
detections are shown in orange, and non-detections are shown in gray. 


et al. (2013) observed the optical transient with the X-Shooter 
instrument on the ESO 8.2-m VLT. In their spectrum ex¬ 
tending from ^3000 to ^ 24000A, they identified several 
weak absorption lines from which they derived a redshift of 
z = 1.874. Both spectra are shown in Figure 4. 

The source was detected by Swift XRT (Page 2013), but 
with insufficient photons for spectral analysis. The source 
was observed with ATCA, but no radio emission was detected. 
Largely because in our sample this is the oldest afterglow at 
the time of discovery, there are not enough broadband data to 
constrain the blast wave physics. 

3.3. GRB131231MPTF13ekl 

GRB 131231A was detected by Fermi LAT (Sonbas et al. 
2013) and GBM (Jenke & Xiong 2014), with photons of en¬ 
ergies up to 9.7 GeV. Xu et al. (2013b) observed the LAT 
error circle with the 1-m telescope at Mt. Nanshan, Xin¬ 
jiang, China. At 7.9 hr after the burst, they detected a single 
R 17.6 mag source that was not present in SDSS images. 
At 17.3 hr after the burst, Malesani et al. (2013) observed the 
afterglow candidate with the MOSaic CAmera (MOSCA) on 
the 2.56-m Nordic Optical Telescope (NOT). The source had 
faded ioR= 18.6. 
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Figure 3. Light curves and SEDs of GBM-iPTF afterglows; GRB 130702A / iPTFlSbxl is shown here in the print version. 
(An extended version of this figure set is available in the online journal.) 
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Table 4 

Log of Spectroscopic Observations 


Date 

Telescope 

Instrument 

Wavelengths (A) 

Lines 

References 


GRB 13101 lA/iPTF13dsw 

2013 Oct 12 08:56 

Gemini South 

GMOS 

5100-9300 

none 

Kasliwal et al. (2013) 

2013 Oct 13 03:59 

ESOA^LT UTS 

X-shooter 

3100-5560 

Lya, Si II, C II, C IV, Al II 

Rau et al. (2013) 




5550-10050 

Fe II, Mg II 


GRB 131231A/iPTF13ekl 

2014 Jan 01 02:15 

Gemini South 

GMOS 

6000-10000 

[OII], [O III], Ca II H+K 

Cucchiara (2014) 

GRB 140508A/iPTF14aue 

2014 May 08 18:55 

HCT 

HFOSC 

3800-8400 

Fe II, Mg II 

Bhalerao & Sahu (2014) 

2014 May 09 06:33 

APO 

DIS 

3200-9800 

none 

none 

GRB 140606B/iPTF14bfu 

2014 Jun 07 19:16 

Keck II 

DEIMOS 

4500-9600 

[OII], [O III], Ua, Ca II H+K 

Perley et al. (2014a) 

GRB 140620A/iPTF14cva 

2014 Jun 20 14:00 

Gemini North 

GMOS 

5090-9300 

Mg I, Mg II, Fe II, Al II, Si II, Si II* 

Kasliwal et al. (2014) 




4000-6600 



GRB 140623A/iPTF14cyb 

2014 Jun 23 08:10 

Gemini North 

GMOS 

4000-6600 

Mg II, Fe II, Al II, Si II, Al III, C I, C IV 

Bhalerao et al. (2014) 

GRB 140808A/iPTF14eag 

2014 Aug 08 21:43 

GTC 

OSIRIS 

3630-7500 

DLA, S II, Si II, O I, C II, Si IV, Fe II, Al II, C IV 

Gorosabel et al. (2014a) 



GS+GMOS (normalized) 

GRB 131231A/iPTF13ekl+0.90 days 
I^ 



z = 0.6419 




(normalized) St I 

lA / iPTF14cvb+0.12 davs I 


GN+GMOS,( 

GRB 140623A / iPTF14cyb +0.12 days ' 



z= 1.06 
z=1.92 




GRB 140808A / iPTF14eag +0.89 days 
fLj _ I \l I _L 

4000 5000 6000 7000 8000 

observed wavelength (A) 


z = 3.293 
J- 


9000 


Figure 4. Afterglow spectra. The horizontal axis shows wavelength in vacuum in the observer frame, and the vertical axis shows scaled flux. Lines at the redshift 
of the putative host are labeled in black; lines corresponding to any intervening absorbing systems are labeled in red. Note that in cases where one or fewer lines 
are discernible in our spectra, the redshifts have been reported in GCNs by other groups. 
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Table 5 

Optical Observations of GBM-iPTF Afterglows 


Date (mid) Inst.^ Mag.^ 


GRB 130702A/iPTF13bxl 


2013 Jul 02 04:18 

P48 

0.18 

R 

17.38± 

0.04 

2013 Jul 02 05:10 

P48 

0.21 

R 

17.52± 

0.04 

2013 Jul 03 04:13 

P60 

1.17 

8 

18.80± 

0.04 

2013 Jul 03 04:15 

P60 

1.17 

i = 

18.42± 

0.04 

2013 Jul 03 06:16 

P60 

1.26 

i = 

18.56± 

0.06 

2013 Jul 03 06:17 

P60 

1.26 

r = 

18.66± 

0.05 

2013 Jul 03 06:20 

P60 

1.26 

8 

18.86± 

0.04 


Note. — A machine readable version of this table is available in the online 
journal. 

^ RATIR data are from Butler et al. (2013b,a, 2014b). GROND data are from 
Sudilovsky et al. (2013). Keck near-infrared data for GRB 140606B are from 
Perley et al. (2014c). 

^ Time in days relative to GBM trigger. 

^ Magnitudes are in the AB system (Oke & Gunn 1983). 


Table 6 

Radio Observations of GBM-iPTF Afterglows 


Date (Start) Inst.^ At^ Flux Density^ 


GRB 130702A/iPTF13bxl 


2013 Jul 04 

CARMA 

2 

±(93) 

_ 

1580 


330 

2013 Jul 04 

VLA 

2.3 

±(5.1) 

= 

1490 


75 

2013 Jul 04 

VLA 

2.3 

±(7.1) 

= 

1600 

± 

81 

2013 Jul 05 

CARMA 

3.1 

±(93) 

= 

1850 

± 

690 

2013 Jul 06 

CARMA 

4.1 

±(93) 

= 

1090 


350 

2013 Jul 08 

CARMA 

6.1 

±(93) 

= 

1440 

± 

260 

2013 Jul 08 

CARMA 

7 

±(93) 

= 

1160 

± 

320 

2013 Jul 14 

CARMA 

12 

±(93) 

= 

900 


230 

2013 Jul 15 

CARMA 

13 

±(93) 

= 

1550 


590 

2013 Jul 24 

CARMA 

22 

±(93) 

= 

1430 

± 

480 

2013 Jul 25 

CARMA 

23 

±(93) 

< 

1890 



2013 Aug 12 

CARMA 

41 

±(93) 

= 

450 


210 


Note. — A machine readable version of this table is available in the online 
journal. 


^ The ATCA observation is from Hancock et al. (2013). 

^ Time in days relative to GBM trigger. 

^ Flux density in /xJy as a function of frequency in GHz. For detections, the 
confidence intervals are Icr statistical uncertainties added in quadrature with 
an estimated 5% systematic error. For non-detections, we show 3a upper 
limits. 

Although we had imaged 10 P48 fields shortly after the 
Fermi trigger (Singer et al. 2013a), due to the short visibility 
window at Palomar we were only able to obtain one epoch. 
At 1.45 hr after the burst, we detected an 7? = 15.7 mag optical 
transient iPTFlSekl at the position of the Nanshan candidate. 
Though our single detection of iPTFlSekl could not by itself 
rule out that the source was a moving solar system object, the 
Nanshan detection at 6.46 hr, fitting a decay with a power-law 
index of a = 1.03, was strong evidence that the transient was 
the optical afterglow of GRB 131231A. 

On January 1.09 UT (21.5 hr after the trigger), we observed 
the afterglow with Gemini South using the GMOS camera 
(Hook et al. 2004) in Nod&Shuffle mode: we obtained 32 
dithered observations of 30 s each at an average airmass of 
2. We analyzed this data set using the dedicated GEMINI 
package under the IRAF environment and extracted the 1- 
dimensional spectrum using the APALL task. We determined 


the redshift of the GRB, based on the simultaneous identi¬ 
fication of forbidden nebular emission lines ([O II], [O III]) 
and absorption features (CaH&K) at the same redshift of 
z = 0.6419. In Figure 4, we show the normalized spectrum. 

The source was also detected by Swift XRT (Mangano et al. 
2014b) and UVOT (Holland & Mangano 2014), as well as 
CARMA (Perley 2014). 

With only the millimeter, optical, and X-ray observations, 
the SED is highly degenerate. Contributing to the degeneracy, 
the X-ray and optical observations appear to fall on the same 
power-law segment. It is consistent with either fast or slow 
cooling if the greater of Uc of is near the optical, assuming 
a fiat electron distribution with p ^ 1.5. It is also consistent 
with slow cooling if Vc is above the X-ray band and p 2.6. 

3.4. GRB 140508A/iPTF14aue 

This burst was detected by Fermi GBM and INTEGRAL 
SPI-ACS (Yu & Goldstein 2014), as well as by Konus-WmJ, 
Mars Odyssey (not included in the GCN circular). Swift 
BAT (outside the coded FOV), and MESSENGER, yielding 
a l.°5 X 12' IPN error box (Hurley et al. 2014c). 

Due to poor weather early in the night, P48 observations 
started 6.7 hr after the trigger (Singer et al. 2014a). We found 
one optical transient candidate within the IPN triangulation, 
iPTF14aue, which faded from r = 17.89 ±0.01 mag with a 
power-law fit of a = 1.12 ± 0.1 over a timescale of 1.5 hr. 

We triggered a Swift TOO. From 0.8 to 8.1 days after 
the trigger. Swift XRT detected a coincident X-ray source 
that faded with a power law a = 1.48 (+0.15, -0.14) (Amaral- 
Rogers 2014a,b). The source was also detected by Swift 
UVOT (Marshall & Amarel-Rogers 2014). 

Moskvitin et al. (2014) obtained a 20-minute, 3800-7200 A 
spectrum of iPTF14aue with the 6-m BTA telescope in Ze- 
lenchukskaia. Exhibiting no absorption features, this estab¬ 
lished an upper limit of z < 2.1. Malesani et al. (2014) 
used the Andalucia Faint Object Spectrograph and Cam¬ 
era (ALFOSC) on NOT to get an 1800 s spectrum spanning 
3200-9100 A, and found several absorption features at red¬ 
shift z = 1.03. Consistent redshifts were reported by Wiersema 
et al. (2014) with the ACAM instrument on the 4.2-m William 
Herschel Telescope and by Bhalerao & Sahu (2014) with Hi¬ 
malaya Faint Object Spectrograph and Camera (HFOSC) on 
the 2-m HCT. This last spectrum is shown in Figure 4. 

Due to the brightness of the optical transient, optical pho¬ 
tometry was available from several facilities up to 4.5 days 
after the burst (Gorosabel et al. 2014b; Moskvitin et al. 2014; 
Malesani et al. 2014; Masi 2014; Butler et al. 2014b,a; Fuji- 
wara et al. 2014; Voinova et al. 2014a). 

Horesh et al. (2014) detected the source with VLA 5.2 days 
after the Fermi trigger, at 6.1 GHz (C band) and at 22 GHz (K 
band). A broadband SED constructed from P60 and XRT data 
from around this time is consistent with p ^2. Because p is 
not distinguishable from 2, we cannot discriminate between 
fast and slow cooling based on this one time slice. However, 
given the late time of this observation, the slow cooling in¬ 
terpretation is more likely, putting Um between the radio and 
optical bands and between the optical and X-ray. Because 
the VLA light curve is decreasing with time, an ISM circum- 
burst density profile is favored. 

3.5. GRB 140606B/iPTF14bfu 

Fermi trigger 423717114 (Burns 2014) was observable 
from Palomar for several hours, starting about 4.3 hr after the 
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time of the burst. Based on the final GBM localization, we 
searched ten P48 fields and found several plausible optical 
transient candidates (Singer et al. 2014c). 

iPTF14bfu had no previous detections in iPTF between 
2013 May 23 and October 13. Its position was outside the 
SDSS survey footprint, but it had no plausible host associa¬ 
tions in VizieR (Ochsenbein et al. 2000). From 4.3 to 5.5 hr 
after the burst, it faded from R = 19.89 ± 0.10 to 20.32 ± 0.14 
mag, fitting a power law of a = -1.6 ± 0.7 relative to the time 
of the GBM trigger. iPTF14bfw (R= 19.96 ± 0.06 mag) was 
coincident with an r = 21.27 galaxy in SDSS DRIO and dis¬ 
played no statistically significant photometric variation over 
the course of our P48 observations. iPTF14bgc (R = 18.44 ± 
0.02 mag) was coincident with an R = 21.07 ± 0.08 mag point 
source in our co-added reference image composed of expo¬ 
sures from 2013 July 31 through September 24. iPTF14bga 
(R= 19.75±0.06 mag) was likewise coincident with a R = 
20.42 ±0.17 mag point source in our reference image com¬ 
posed of exposures from 2011 July 29 through October 20. 

On the following night, we observed all four candidates 
again with P48 and P60 (Perley & Singer 2014). iPTF14bfw 
and iPTF14bga had not faded relative to the previous night. 
iPTF14bgc had faded to R = 20.68 ± 0.21 mag, consistent 
with the counterpart in our reference images but significantly 
fainter than the previous night. A power-law fit to the decay 
gave a temporal index of a = -l.l±0.1, entirely consistent 
with typical GRB afterglows. iPTF14bfu was not detected in 
our P48 images to a limiting magnitude of R < 21.1, but it 
was detected in stacked P60 images (r = 21.1 ± 0.2), consis¬ 
tent with a power law of a ^ -0.5. 

An IPN triangulation from Fermi, Konus-JFmJ, and MES¬ 
SENGER yielded a long, slender 14.° 18 x 0.°414 error box 
that contained iPTF14bfu and iPTF14bfw (Hurley et al. 
2014d). 

We obtained two 900 s spectra with the DEIMOS spectro¬ 
graph on the Keck II 10 m telescope (Perley et al. 2014a). On 
a blue continuum, we found [O II], [O III], and H a emission 
features and Ca II absorption features, at a common redshift 
of z = 0.384. A galaxy offset by ^ 2" along the slit showed 
the same emission lines at the same redshift. 

Swift XRT observed the location of iPTF14bfu for a total of 
9 ks from 2.1 to 9.3 days after the GBM trigger, and found a 
source that faded with a power-law fit of a = -1.0 (-1-0.7, -0.6) 
(Mangano et al. 2014a; Mangano & Burrows 2014; Mangano 
2014). 

At 18.4 days after the trigger, we obtained a 1200 s spec¬ 
trum of iPTF14bfu with the Low Resolution Imaging Spec¬ 
trometer (LRIS) on the Keck I 10-meter telescope (Perley 
et al. 2014b). The spectrum had developed broad emission 
features. A comparison using Superfit (Howell et al. 2005) 
showed a good match to SN 1998bw near maximum light, in¬ 
dicating that the source had evolved into an SNIc-BL. Our 
late-time photometry and spectroscopy will published sepa¬ 
rately (Cano et al. 2015). 

Although there were three radio detections of 
GRB 140606B, only during the first CARMA detection 
does the optical emission appear to be dominated by the 
afterglow. We can construct an SED around this time using 
nearly coeval DCT and XRT data. Because of the faintness 
of the X-ray afterglow, the spectral slopes /3x and Pox are 
only weakly determined. As a result, there is a degeneracy 
between two plausible fits. The first has Um anywhere below 
the CARMA band, Vc just below the X-rays, and p ^2. The 
second has Um just above the radio and Uc in the middle of the 


XRT band, with p ^2.2. 

The early P48 observations do not connect smoothly with 
the P60 and DCT observations from Ar = 1 to 4 days. This 
may indicate a steep-shallow-steep time evolution requiring 
late-time energy injection, or may just indicate that the after¬ 
glow is contaminated by light from the host galaxy or the SN 
at relatively early times. 

3.6. GRB 140620A/iPTF14cva 

This burst is distinctive in our sample for two reasons. First, 
it is the earliest afterglow detection in the iPTF sample at /St = 
0.25 hr. Second, its broadband SED is not clearly explainable 
by the standard forward shock model. 

Fermi trigger 424934131 (Eitzpatrick & Connaughton 
2014) was observable from Palomar for about 6 hr from the 
time of the burst. Based on the ground localization, we started 
observing ten P48 fields about 10 minutes after the trigger. 
Based on the final localization, we added 10 more fields, for 
a total of 20, about an hour after the trigger. 

The candidate iPTE14cva (Kasliwal et al. 2014) was con¬ 
tained within one of the early 10 fields. Erom 14.9 to 87.2 
minutes after the trigger, the candidate faded from R = 17.60 ± 
0.01 to 18.80 ±0.02 mag, consistent with a somewhat slow 
power law of a = 0.62 ±0.01. 

We observed the candidate with CMOS on the 8-m Gemini 
North telescope. Starting 8.8 hr after the trigger, we obtained 
two 900 s spectra extending from 4000 to 9300 A. We de¬ 
tected Mg II and Ee II absorption lines at z = 0.88 and many 
absorption features at a common redshift of z = 2.04. The 
lack of Fya absorption implied an upper limit of z ^ 2.3 and 
suggested that z = 2.04 was the redshift of the source. 

We triggered Swift and VLA follow-up. In a 3 ks exposure 
starting 10.4 hr after the Fermi trigger. Swift XRT detected an 
X-ray source with a count rate of 1.2 x 10“^ counts s“^ (De 
Pasquale 2014b). Over the next four days of Swift observa¬ 
tions, the X-ray source faded with a slope a= 1.32 ±0.16 
(De Pasquale 2014a). A fading source was also detected by 
Swift UVOT (Siegel & De Pasquale 2014). 

The source was detected by VLA on June 23 at 6.1 GHz (C 
band) at 108 ±15 /iJy and at 22 GHz {K band) at 62 ± 15 /iJy. 
On June 30, there was a marginal detection in C band with 
48 ± 12 /iJy and no detection in K band with a noise level of 
15 /iJy rms. 

The optical transient was also observed in R band by the 
Konkoly Observatory (Kelemen 2014) and the 1-m telescope 
at the Tien Shan Astronomical Observatory (Voinova et al. 
2014b). 

The SED of this afterglow cannot be explained by a stan¬ 
dard forward shock model. If we place the peak frequency 
near the radio band, the optical and X-ray fluxes are drasti¬ 
cally underpredicted, whereas if we place the peak frequency 
between the optical and X-ray bands, we miss the radio ob¬ 
servations by orders of magnitude. This seems to require an 
additional component. One possibility is that there is a for¬ 
ward shock peak in the UV and a reverse shock peak at low 
frequencies (similar to GRB 130427A; see Laskar et al. 2013; 
Perley et al. 2014d). Another possibility is that there is an 
inverse Compton peak in the UV (similar to GRB 120326A; 
Urata et al. 2014). 

3.7. GRB 140623A/iPTFMcyb 

Fermi trigger 425193729 (von Kienlin 2014) was observ¬ 
able from Palomar for about 6 hr from the time of the burst. 


The Needle in the 100 deg^ Haystack 


13 


Based on the ground localization, we started imaging 10 fields 
11 minutes after the trigger. The final Fermi localization, 
which was avilable 2.6 hr later, shifted by 13.°4. Due to the 
large change in the localization, we calculated only a 4% 
chance that the source was contained within the P48 fields. 

Candidate iPTF14cyb (Kasliwal et al. 2014), situated at an 
extreme edge of the P48 tiling, was within the Icr confidence 
region for both the ground and final localizations. From 16 to 
83 minutes after the trigger, the source faded from /?= 18.04 ± 
0.01 to 19.69 ± 0.06 mag, consistent with a power-law decay 
with an index a = 0.94 ± 0.03. 

Starting 2.8 hr after the trigger, we obtained two 900 s 
CMOS spectra extending from 4000 to 9300 A. We detected 
Mg II and Fe II absorption lines at z = 1.06 and many absorp¬ 
tion features at z = 1.92. The lack of Fya absorption implied 
that this was the redshift of the burst. 

We triggered Swift, VLA, and CARMA follow-up. In a 3 ks 
exposure starting 10.7 hr after the burst. Swift XRT detected 
an uncataloged X-ray source with a count rate of (2.2 ± 0.6) x 
10“^ counts s“^ (D’Elia et al. 2014). By 79 hr after the trigger, 
the source was no longer detected in a 5 ks exposure (D’Elia 
& Izzo 2014). No radio source was detected with VLA in 
C band (6.1 GHz) to an rms level of 17 /iJy, or in K band 
(22 GHz) to an rms level of 18 /iJy. 

Because of the lack of radio detections and the extreme 
faintness of the X-ray afterglow, the broadband behavior of 
the afterglow does not constrain the shock physics. 

3.8. GRB 140808A/iPTF14eag 

Fermi trigger 429152043 (Zhang 2014) was observable 
from Palomar about 3 hr after the burst. We imaged 13 
fields with P48 and found one promising optical transient. 
iPTE14eag was situated on the extreme edge of one of the P48 
tiles that was just outside the GBM la contour. It faded from 
R = 18.91 ±0.06 to 19.29 ±0.10 mag from 3.35 to 4.91 hr af¬ 
ter the trigger and had no archival counterparts in SDSS or in 
our own reference images. 

Due to Hurricane Iselle, we were unable to use our TOO 
programs on Keck or Gemini North. We requested photomet¬ 
ric confirmation of the fading from HCT (Sahu et al. 2014), 
submitted a Swift TOO, and sent our GCN circular (Singer 
et al. 2014b) to encourage others to obtain a spectrum. 

Swift observed the position of iPTE14eag from 
11.6 to 14.4 hr after the burst (Page et al. 2014). An X-ray 
source was detected with a count rate of 1.5 x 10“^ counts s“^. 
In a second observation starting 62.2 hr after the trigger 
(Page & Cenko 2014), the source had faded to below 
2.46 X 10“^ counts s“^. No source was detected by UVOT 
(Oates & Cenko 2014). 

We obtained a spectrum with the OSIRIS instrument (Cepa 
et al. 2000) on Gran Telescopio Canarias (GTC) in 4 x 900 s 
exposures with a mean epoch of 21.340 hr after the burst. We 
used the RIOOOB grism and a 1" slit, with a resolution of 
R 1000. We determined a redshift of 3.293 (improved from 
the one given by Gorosabel et al. 2014a) through the iden¬ 
tification of strong absorption features. The fiux-calibrated 
spectrum is shown in Eigure 4. 

The source was detected in radio with VLA (Corsi & 
Horesh 2014) and AMI (Anderson et al. 2014a). The broad¬ 
band SED around the time of the VLA detection broadly fits 
a forward shock model but is poorly constrained due to the 
lack of a contemporaneous X-ray detection. The spectral 
slope between the two VLA bands is somewhat steeper than 


the standard low-frequency value of /3 = -1/3, possibly indi¬ 
cating that the radio emission is self-absorbed. We obtained 
14 AMI observations every 2 or 3 days from 2014 August 8 
until 2014 September 12. Observations were 2-4 hr in du¬ 
ration. AMI first detected the afterglow 4.6 days post-burst. 
The AMI light curve peaked ^10.6 days post-burst at 15.7 
GHz, which is characteristic of forward shock emission at ra¬ 
dio wavelengths (Chandra & Erail 2012). 

A peculiar feature of the optical light curve is that the P60 r- 
and /-band observations at Ar 2 days appears to be inverted, 
with a rising rather than falling spectral shape, compared to 
the earlier P60 photometry at Ar 1 day. However, this fea¬ 
ture is within the error bars and may be merely a statistical 
fiuctuation. 

This is the highest-redshift burst in our sample and also had 
the weakest prompt emission in terms of the fiuence observed 
by GBM. 

4. THE POPULATION IN CONTEXT 
4.1. Selection Effects 

Eirst, we investigate the properties of the subset of GBM 
bursts followed up by iPTE compared to the GBM bursts as a 
whole. It is known that, on average, GRBs with larger prompt 
fiuences have brighter optical afterglows, though the corre¬ 
lation is very weak (Nysewander et al. 2009). In Eigure 6, 
we plot the fiuence in the 10-1000 keV band and 1-a local¬ 
ization radius of all GBM bursts from the beginning of our 
experiment, retrieved from the Fermi GBM Burst Catalog at 
HEASARC."^^ As expected, there is a weak but clearly dis¬ 
cernible correlation between fiuence and radius, F oc 
with a Pearson correlation coefficient of R = 0.64."^^ The sub¬ 
set of bursts that we followed up spans a wide range in fiuence 
and error radii up to ^ 10°. The bursts for which we detected 
optical afterglows are preferentially brighter, with the faintest 
burst having a fiuence as low as 3 x 10“^ ergcm“^. There 
are some bright (> 3 x 10“^ ergcm“^) and well-confined 
(< 1.8°) events for which we did not find afterglows: those 
at 2013 August 28 07:19:56, 2013 November 27 14:12:14, 
and 2014 January 04 17:32:00 (see Table 3). However, 
these non-detections are not constraining given their ages 
of 20.28, 13.46, and 18.57 hours respectively. Conversely, 
there were two especially young bursts (followed up at /St = 
0.11 and 0.17 hr) for which we did not detect afterglows. 
The non-detection of the burst at 2014 July 16 07:20:13 
makes sense because we searched only 28% of the GBM 
localization. The non-detection on 2014 Paril 04 04:06:48, 
for which we observed 69% of the localization, is a little 
more surprising, especially given its relatively high fiuence 
of 8 X 10“^ ergcm“^; this is a possible candidate for a “dark 
GRB.” On the whole, however, we can see that (1) we have 
followed up bursts with a large range of error radii and fiu¬ 
ences, (2), there is a weak preference toward detecting bursts 
with small error radii, and (3) the detections tend toward 
bursts with high fiuences. Naively one might expect higher 
fiuences to translate into lower redshifts, but the interplay be¬ 
tween the GRB luminosity function and detector threshold 
greatly complicates such inferences (Butler et al. 2010). 

http://heasarc.gsfc.nasa.govAV3Browse/fermi/ferniigbrst.html 

In a separate sample of GBM GRBs compiled by Connaughton et al. 
(2015), the correlation between error radius and photon fiuence is slightly 
stronger than the correlation between error radius and fiuence. However, we 
use fiuence rather than photon fiuence here because the latter is not available 
for all bursts in the online Fermi GBM archive. 
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GRB 130702A/iPTF13bxl 



GRB 131231A/iPTF13ekl 



co-added reference 




co-added reference 





GRB 140606B / iPTF14bfu co-added reference 









co-added reference 




Figure 5. Gamma-ray localizations, P48 tiles, and discovery images for the GBM-iPTF afterglows. The Fermi GBM 1- and 2-(7 regions are shown as black 
contour lines, the P48 tiles as gray rectangles, the 3cr IPN triangulations in blue (when available), and the LAT Icr error circles in green (when available). The 
positions of the optical transients are marked with black diamonds. 
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GBM radius (°) 

Figure 6. Fluence and statistical error radius of GBM bursts. Orange dots 
mark bursts that were followed up with iPTF; black circles around orange 
dots mark bursts whose afterglows were detected by iPTF. The black line is 
a power-law ht. 



redshift 


Figure 7. Cumulative distribution of redshifts of long GRBs observed by 
Swift BAT (gray) and the GBM-iPTF experiment (orange). 

Second, the rich sample of all of the GRB afterglows that 
we have today is undeniably the result of the success of the 
Swift mission. It is therefore interesting to consider how the 
GBM-iPTF sample is similar to or different from the Swift 
sample, given the differences in bandpasses and our increased 
reliance on the optical afterglow. In Figure 7, we plot the 
cumulative redshift distribution of our sample, alongside the 
distribution of redshifts of long GRBs detected by Swift In¬ 
deed, we find that our sample is at lower redshifts; the former 
distribution lies almost entirely to the left of the latter, and 
the ratio of the median redshifts (z = 1.5 versus z = 1.9) of the 
two populations is about 0.75. However, with the small sam¬ 
ple size, the difference between the two redshift distributions 
is not significant: a two-sample Kolmogorov-Smimov test 
yields a /7-value of 0.26, meaning that there is a 26% chance 
of obtaining these two empirical samples from the same un¬ 
derlying distribution. More GBM-iPTF events are needed to 
determine whether the redshift distribution is significantly dif¬ 
ferent. 

This sample was extracted from the Swift GRB Table, http://swift.gsfc. 
nasa.gov/archive/grb_table/. 


4.2. GRBs as Standard Candles? 

Amati et al. (2002) pointed out a striking empirical correla¬ 
tion in the rest-frame prompt emission spectra of BeppoSAX 
GRBs, with the peak energy (in the uFy sense) ^'peak related 
to the bolometric, isotropic-equivalent energy release £’iso by 
£’peak oc £’iso^- H was quickly realized that such a relation, if 
intrinsic to the bursts, could be used to measure the redshifts 
of GRBs non-spectroscopically (Atteia 2003). As with the 
Phillips relation for SNe la (Phillips 1993), with such a rela¬ 
tion GRBs could serve as standardizable candles in order to 
measure cosmological parameters (Dai et al. 2004; Friedman 
& Bloom 2005; Ghirlanda et al. 2006; etc.). 

However, there has been a vigorous debate about whether 
the Amati relation and related correlations are innate to GRBs 
or reflect a detector-dependent selection bias (Band & Preece 
2005; Ghirlanda et al. 2005; Nakar & Piran 2005; Sakamoto 
et al. 2006; Butler et al. 2007; Cabrera et al. 2007; Schaefer & 
Collazzi 2007; Butler et al. 2009; Firmani et al. 2009; Krimm 
et al. 2009; Butler et al. 2010; Shahmoradi & Nemiroff 2011; 
Collazzi et al. 2012; Kocevski 2012). One alternative inter¬ 
pretation is that bursts to the upper-left boundary of the Amati 
relation are selected against by photon-counting instruments 
because, being relatively hard, there are fewer photons. The 
lack of bursts to the lower right of the Amati line may be due 
to a genuine lack of relativistic explosions that are much softer 
than, but as energetic as, standard GRBs. 

It has been difficult to directly test the Amati relation in 
the context of Fermi bursts because most lack known red¬ 
shifts, since bursts that were coincidentally observed and lo¬ 
calized by the Swift BAT do not directly sample the selec¬ 
tion bias of Fermi GBM. However, Heussaff et al. (2013) 
showed that many Fermi bursts that lack known redshifts 
would be inconsistent with the Amati relation at any dis¬ 
tance. (See also Urata et al. 2012 for outlier events detected 
by Fermi LAT and Suzaku WAM.) Here, we have a small 
sample of Fermi bursts with known redshifts. One of them, 
GRB 140606B/iPTF14bfu at z = 0.384, is a clear outlier, over 
2(7 away from the mean Amati relation. This burst is not 
alone: in Figure 8, we have marked a selection of previ¬ 
ous long GRBs with spectroscopically identified SNe. Three 
among them are also outliers. (A possible caveat is that the 
prompt emission mechanism for GRB 140606B could be dif¬ 
ferent from typical cosmological bursts; we explore this in the 
next section.) To be sure, most of the bursts in our GBM-iPTF 
sample fall within a la band of the Amati relation. This in¬ 
cludes the nearest event to date, GRB 130702A/iPTF13bxl at 
z = 0.145. However, the one outlier in our admittedly small 
sample strengthens the case that the boundary of the Amati 
relation is somewhat influenced by the detector thresholds and 
bandpasses. 

4.3. Shock Breakout 

Two GRBs in our sample, GRB 130702A/iPTF13bxl and 
GRB 140606B/iPTF14bfu, have £'iso ^10^^ erg (rest frame), 
energetically intermediate between “standard” luminous, cos- 
mically distant bursts and nearby llGRBs. Prototypes of the 
latter class include GRB980425/SN 1998bw (Galama et al. 
1998; Kulkami et al. 1998), which was also the first SN dis¬ 
covered in association with a GRB. They offer an interest¬ 
ing test case for competing theories to explain the wide range 
of prompt gamma-ray energy releases observed from GRBs 
(e.g., Schulze et al. 2014b). 

It has been suggested that the two luminosity regimes cor- 
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Figure 8. Rest-frame energetics of GBM-iPTF bursts (in orange) in comparison to an illustrative sample of previous GRB-SNe (in black; includes 
GRB 060218/SN2006aj, Plan et al. 2006, Modjaz et al. 2006, Sollerman et al. 2006; GRB 100316D/SN2010bli, Chornock et al. 2010, Bufano et al. 2011; 
GRB 120422A/SN2012bz, Melandri et al. 2012, Schulze et al. 2014a; GRB 130215A/SN2013ez, Cano et al. 2014; GRB 130427A/SN2013cq, Xu et al. 2013d; 
and GRB 130831A/SN2013fu, Klose et al. 2013). A general long GRB sample from Amati (2006) and Amati et al. (2008, 2009) is shown in gray. The solid 
black line represents the Amati relation as given in Amati (2006), £’peak = 95(£^-),pso/10^^ erg)^-^^ keV. The black dashed lines show the relation’s Icr dispersion 
of ±0.4 dex. 


respond to different prompt emission mechanisms (Bromberg 
et al. 2011). The llGRBs could be explained by the break¬ 
out of a mildly relativistic shock from the progenitor enve¬ 
lope (Nakar & Sari 2012). High-luminosity bursts, on the 
other hand, are thought to be produced by internal shocks 
within an ultra-relativistic jet (Rees & Meszaros 1994) that 
has successfully punched through the star. A central engine 
that sometimes fails to launch an ultra-relativistic jet is one 
way to unify the luminosity functions of standard GRBs and 
llGRBs (Pescalli et al. 2015). 

The smoking gun for the relativistic shock breakout model 
is a cooling, thermal component to the prompt X-ray 
emission, as in the case of GRB 060218 (Campana et al. 
2006). Unfortunately, this diagnostic is not possible for 
GRB 130702A and GRB 140606B because we lack early-time 
Swift observations. 

However, Nakar & Sari (2012) propose a closure relation 
(their Equation (18)) between the prompt energy, temperature, 
and timescale that is valid for shock breakout-powered GRBs. 
We reproduce it here: 


y:0bs 
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(3) 


If we very crudely assume that all of the prompt emission is 
from a shock escaping from the progenitor envelope, then we 
can use E’iso, E’peak, and T^o as proxies for those observables. 
This gives us a simple discriminator of which bursts are plau¬ 
sible shock breakout candidates, the ratio 

^ = (1+z)C/7’90, (4) 

which should be close to 1. As expected, most of the ener¬ 
getic (E’iso > 10^^ erg), cosmic (z > 0.5) GRBs in our sam¬ 
ple are inconsistent with the closure relation. They are all 
much shorter in duration, given their 7 -ray spectra, than 
would be expected for a shock breakout. The exception is 
GRB 140623A/iPTF14cyb, which yields ^ = 0.5 ± 0.5. In 
this particular case, one possible explanation is that the cen¬ 
tral engine simply remained active for much longer than the 
timescale of the shock breakout. 

Surprisingly, of the two low-luminosity, low-redshift bursts 
in our sample, GRB 130702A/iPTF13bxrs prompt emission 
was much too brief to be consistent with this shock break¬ 
out model, with ^ = (1.6 ±0.7) x 10^. Most likely, this 
means that the prompt emission of GRB 130702A is simply 
a very soft, very sub-luminous version of an otherwise ‘ordi¬ 
nary’ long GRB. Any early-time shock breakout signature, if 
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present, was not observed either because it occurred at ener¬ 
gies below GBM’s bandpass or because it was much weaker 
than the emission from the standard GRB mechanism. How¬ 
ever, GRB 140606B/iPTF14bfu’s prompt emission is consis¬ 
tent with the closure relation, with ^ = 0.5 ± 0.3. Though we 
must interpret this with caution because we cannot disentan¬ 
gle a thermal component from the GBM data, if we naively 
apply linear least squares to (the logarithm of) Equations (14), 
(16), (17) of Nakar & Sari (2012), 


Eho - 

s 2 X erg, 

(5) 

Tbo ~ 507/^0 keV, 

(6) 

^obs ^ 
±0 

e 10^ S, 

(7) 


V.o 



then we find the breakout radius and Lorentz factor to be 

7?bo = (l-3±0.2)x 10^7?©, 

7 /,o = 14±2. 

The breakout radius is comparable to that which Nakar & 
Sari (2012) find for GRB 060218 and GRB 100316D, sug¬ 
gestive of breakout from a dense wind environment, rather 
than the star itself.^^ However, the derived Lorentz factor of 
GRB 140606B is a bit higher than those of the other two ex¬ 
amples. 

Another way to constrain the nature of the explosion is to 
look at the kinetic energy iso of the blast compared to the 
promptly radiated energy £’^,iso = £’iso and the radiative effi¬ 
ciency T] = £’^^iso/(£’jb,iso +^iso)- After the end of any plateau 
phase, the X-ray fiux is a fairly clean diagnostic of 
assuming that the X-rays are above the cooling frequency 
(Freedman & Waxman 2001). During the slow-cooling phase 
and under the typical conditions where p ^ 2 and 
the X-ray fiux is only weakly sensitive to global parameters 
such as the fraction of the internal energy partitioned to elec¬ 
trons and to the magnetic (e^, e^). Even the radiative losses, 
necessary for extrapolating from the late-time afterglow to the 
end of the prompt phase, are minor, amounting to order unity 
aX At = 1 day (Lloyd-Ronning & Zhang 2004). We calculate 
the isotropic-equivalent rest-frame X-ray luminosity from the 
fiux at Ar = 1 day using Equation (1) of Racusin et al. (201 1), 
reproduced below: 

Lxit) = 47rZ)L Vx(r)(l +2 )-«x+/3x-1 _ (8) 

Then we estimate the kinetic energy at the end of the prompt 
emission phase using Equation (7) of Lloyd-Ronning & 
Zhang (2004): 
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The correction factor R for radiative losses is given by Equa- 


Note that SN 2008D, which seems to be the only case so far of shock 
breakout observed in an “ordinary” SN Ibc, had a 500 s emission episode 
that was not strictly consistent with the picture of shock breakout from a 
progenitor envelope. Svirski & Nakar (2014a, b) explore the case of shock 
breakout through a thick Wolf-Rayet wind, which can accommodate longer 
emission timescales. 


tion (8) of Lloyd-Ronning & Zhang (2004), adopted here: 


R = 


k) 


(17/16)e, 


( 10 ) 


The numeric subscripts follow the usual convention for rep¬ 
resenting quantities in powers of 10 times the cgs unit, i.e., 
e,-i = €b -2 = eb/IO-^, and I'ls = r//(10‘* Hz). We 

assume = 0.1 and cp = 0.01. Eor bursts that have XRT 
detections around At = I day (GRBs 130702A, 131231A, 
140508A, 140606B, and 140620A), we calculate Lx by in¬ 
terpolating a least-squares power-law fit to the X-ray light 
curve. Some of our bursts (GRBs 131011 A, 140623A, and 
140808A) were only weakly detected by XRT; for these we 
extrapolate from the mean time of the XRT detection assum¬ 
ing a typical temporal slope of ax = L43 ±0.35 (Racusin 
et al. 2011). The kinetic and radiative energies of our eight 
bursts are shown in Eigure 9. Half of our bursts are reason¬ 
ably well constrained in E^-E^ space; these are shown as red 
points. The other half (GRBs 131011 A, 131231A, 140620A, 
and 140623A) have highly degenerate SEDs, so their position 
in this plot is highly sensitive to model assumptions; these 
are shown as gray points. Dotted lines are lines of constant 
radiative efficiency. 

Within our sample, there are at least three orders of magni¬ 
tude of variation in both Ek^i^o and £’^,iso- The two GRB-SNe 
have radiative and kinetic energies of ^ 10^^ erg, both two to 
three orders of magnitude lower than the other extreme in our 
sample or the average values for Swift bursts. In our sample, 
they have two of the lowest inferred radiative efficiencies of 
7 ^ ^ 0.1-0.5, but these values are not atypical of BATSE bursts 
(e.g., Lloyd-Ronning & Zhang 2004) and are close to the me¬ 
dian value for Swift bursts. These are, therefore, truly less 
energetic than cosmological bursts, not merely less efficient 
at producing gamma rays. 


5. LOOKING EORWARD 

In this experiment, we have followed up 35 Eermi GBM 
bursts, scanning areas from 30 to 147 deg^. To date, we have 
detected eight afterglows with apparent optical magnitudes as 
bright as R 16 and as faint as R ^ 20. We have found red- 
shifts as nearby as z = 0.145 and as distant as z = 3.29. A 
continuation of the project should reveal more low-redshift 
events, more GRB-SNe, and more relatively hard GRBs. 

We aim to uncover the much fainter afterglows of short, 
hard bursts by using stacked P48 exposures and integrating a 
co-addition stage into the real-time pipeline, and by honing 
our follow-up to sift through the increased number of candi¬ 
dates. The greatest factor limiting discoveries is, of course, 
that Eermi detects bursts all over the sky, only a fraction of 
which are visible from Palomar. Given our success so far, 
we enthusiastically suggest that other wide-field surveys im¬ 
plement a similar program. Eurthermore, automatically shar¬ 
ing lists of candidates between longitudinally separated in¬ 
struments would facilitate rapid identification and follow-up 
of the fastest-fading events. 

It is uncertain what directions future gamma-ray space mis¬ 
sions will take. Some may be like Swift, able to rapidly 
train multiple on-board follow-up instruments on new targets. 
Even if they lack these capabilities, we should be able to rou¬ 
tinely locate GRB afterglows and find their redshifts using tar¬ 
geted, ground-based optical transient searches similar to the 
one that we have described. 
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will present no undue difficulties for ZTF when we consider 
its 15-fold increase in areal survey rate as compared to iPTF. 

However, there is a second challenge that these optical sig¬ 
natures are predicted to be fainter than perhaps 22 mag (with 
the exception of on-axis afterglows, which should be rare but 
bright due to beaming). For meter-size telescopes, this will 
require integrating for much longer (10 minutes to 1 hr) than 
we have been performing with iPTF. Fortunately, because the 
LIGO antenna pattern is preferentially sensitive above and di¬ 
rectly opposite of North America, we are optimistic that many 
early Advanced LIGO events should be promptly accessible 
from Palomar with long observability windows (Kasliwal & 
Nissanke 2014). 

The main difficulty for any GW optical counterpart search 
will be the inundation of false positives due to the required 
depth and area. We enumerate the following strategies to help 
identify the one needle in the haystack: 

1. Improved machine learning algorithms (see Bue et al. 
2014 in the context of iPTF) will decrease the contam¬ 
ination of the discovery stream by artifacts. 


^'y,iso 

Figure 9. Fireball kinetic energy E^^iso at r = Tgo as estimated from X-ray 
flux vs. rest-frame isotropic-equivalent gamma-ray energy Red points 

denote bursts for which E^^iso can be reliably estimated from the Swift XRT 
data; gray points denote bursts for which the calculation of E^^iso naay have 
extreme model dependence. Dashed lines are lines of constant radiative effi¬ 
ciency r] = Ej^[so/(Ek,iso +^7,iso)- The gray, blue, and red rectangles show the 
Icr parameter ranges of Swift BAT, BAT+GBM, and BAT+LAT long GRBs 
from Racusin et al. (2011). 


Looking beyond GRBs, our present effort serves as a proto¬ 
type for searching for optical counterparts of GW transients. 
We expect that many of the techniques that we have described 
and the lessons that we have learned in the context of iPTF 
will generalize to other wide-field instruments on meter-class 
and larger telescopes. 

Near the end of 2015, Advanced LIGO will begin taking 
data, with Advanced Virgo soon following suit. The first bi¬ 
nary neutron star merger detections are anticipated by 2016 
or later (Aasi et al. 2013). On a similar timescale, iPTF will 
transform into the Zwicky Transient Facility, featuring a new 
47 deg^ survey camera that can reach R = 20.4 mag in 30 s. 
The prime GW sources, binary neutron star (BNS) merg¬ 
ers, may also produce a variety of optical transients: on- or 
off-axis afterglows (van Eerten & MacFadyen 2011; Urata 
et al. 2015), kilonovae (Li & Paczynski 1998; Barnes & Kasen 
2013), and neutron-powered precursors (Metzger et al. 2015); 
see Figure 10 for some examples. 

There will be two key challenges. First, GW localizations 
can be even coarser than Fermi GBM error circles. Starting 
around ^600 deg^ in the initial (2015) two-detector configu¬ 
ration (Kasliwal & Nissanke 2014; Singer et al. 2014e), the 
areas will shrink to ^200 deg^ with the addition of Virgo in 
2016. They should reach ^10 deg^ toward the end of the 
decade as the three detectors approach final design sensitivity 
and can approach ~ 1 deg^ as additional planned GW facili¬ 
ties come online (LIGO-India and KAGRA; see Schutz 2011; 
Veitch et al. 2012; Fairhurst 2014; Nissanke et al. 2013; Aasi 
et al. 2013). Since the detection efficiency of our GBM-iPTF 
afterglow search is consistent with the areas that we searched, 
we expect that even the earliest Advanced LIGO localizations 


2. Combining a catalog of nearby galaxies with the dis¬ 
tance and position information from the GW observa¬ 
tions can help to reduce and prioritize targets for further 
follow-up (Nissanke et al. 2013). 

3. Better leveraging of light-curve history across multiple 
surveys will help to automate the selection of targets for 
photometric follow-up with multiple telescopes. 

4. Our first experiences with detections and non¬ 
detections will guide decisions about the optimal filter. 
At the moment, kilonova models prefer redder filters 
(suggesting i band), and precursor models prefer bluer 
(suggesting g band). 

The combination of gamma-ray missions, ground-based 
GW detectors, and synoptic optical survey instruments is 
poised to make major discoveries over the next few years, of 
which we have provided a small taste in this work. We offer 
both lessons learned and a way forward in this multimessen¬ 
ger effort. The ultimate reward will be joint observations of a 
compact binary merger in gamma. X-rays, optical, and GWs, 
giving us an exceptionally complete record of a complex as- 
trophysical process: it will be almost as good as being there. 
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Figure 10. Light curves of short GRB afterglows, scaled to an Advanced 
LIGO range of 200 Mpc. Thin gray lines are afterglows of Swift short GRBs 
that have known redshifts. Thick black lines are synthetic on-axis afterglows 
from van Eerten & MacFadyen (2011) with jet half-opening angles of 0.2 rad 
and observer angles of 0 rad. Jet energies Ej, in units of erg, and circumburst 
densities n, in units of cm“^, are labeled on the plot. The solid, deep-red line 
is the r-band neutron-powered kilonova precursor model from Metzger et al. 
(2015) with opacity Kr = 30 cm^ g“k free neutron mass = 10“^ Mq, and 
electron fraction Ye = 0.05. The solid, light-red lines represent kilonova mod¬ 
els from Barnes & Kasen (2013) with ejected masses of mej = 10“^, 10“^, or 
10“^ Mq and characteristic velocities /3 = v/c = 0.1, 0.2, or 0.3. The Metzger 
et al. (2015) kilonova precursors are blue: they peak at about 0.1 mag brighter 
in the g band than in r. The Barnes & Kasen (2013) kilonova models are red: 
they are about 1 mag brighter in the i band than in r. 
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APPENDIX 

We illustrate three stages of the iPTF pipeline that we dis¬ 
cussed in Section 2: the TOO Marshal (Figure 11), the Trea¬ 
sures Portal (Figure 12), and the Transient Marshal (Fig¬ 
ure 13). 
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Figure 11. Screenshot of the iPTF TOO Marshal shortly after a Fermi GBM 
detection. At this stage, the application presents the recommended P48 fields, 
the time window of observability, and the history of GCN notices and circu¬ 
lars related to the trigger. It gives the human participants the option to cus¬ 
tomize the P48 sequence by adding or removing P48 fields and tuning the 
airmass limit, cadence, or number of images. 
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Figure 12. Screenshot of the Treasures portal, showing new, reference, sub¬ 
traction, and archival SDSS images as well as P48 light curves. This page is 
for the date and field containing GRB 130702A / iPTFlSbxl. 
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Figure 13. Screenshot of the iPTF Transient Marshal, showing 
GRB 130702A/iPTF13bxl. 
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